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Abstract. We investigate a pair creation cascade in the 
magnctosphcre of a rapidly rotating neutron star. We solve 
the set of the Poisson equation for the electro-static poten- 
tial and the Boltzmann equations for electrons, positrons, 
and gamma-ray photons simultaneously. In this paper, we 
first examine the time-dependent nature of particle ac- 
celerators by solving the non-stationary Boltzmann equa- 
tions on the two-dimensional poloidal plane in which both 
the rotational and magnetic axes reside. Evaluating the 
temperature of the heated polar cap surface, which is lo- 
cated near the magnetic pole, by the bombardment of gap- 
accelerated particles, and applying the scheme to millisec- 
ond pulsar parameters, we demonstrate that the solution 
converges to a stationary solution of which pair-creation 
cascade is maintained by the heated polar-cap emission, 
in a wide range of three-dimensional parameter space (pe- 
riod, period derivative, magnetic inclination angle). We 
also present the deathlines of millisecond pulsars. 



1. Introduction 

The Energetic Gamma-Ray Experiment Telescope 
(EGRET) on the Compton Gamma Ray Observatory 
(CGRO) had detected gamma-ray pulsed emissions 
from seven pulsars. The light curves and spectra of 
these gamma-ray pulsars enable us to understand the 
7-ray emission mechanism and examine different models. 
The prospective 7-ray telescope such as Gamma-ray 
Large Area Space Telescope (GLAST) will provide more 
7-ray observation with higher sensitivity. The origin of 
gamma-ray emission from pulsar magnetosphere will be 
studied more thoroughly. 

Base on the work Cheng, Ho, & Ruderman (1986a, b), 
the three-dimensional outer gap models (Romani & Yadi- 
garoglu 1995; Cheng, Ruderman and Zhang 2000) can suc- 
cessfully explain the observed light curve in reproducing 
the wide separation of the two peaks with bridge emis- 
sion. Dyks and Rudak (2003) found that radiation region 
with inner boundary extended toward the stellar surface 
can explain the off-pulse emission, which cannot be re- 
produced by a traditional outer gap. An one-dimensional 



model of an outer gap with electrodynamical approach 
was investigated (Hirotani & Shibata 1999a, b, c; Hirotani, 
Harding & Shibata 2003). They solved accelerating elec- 
tric field, curvature radiation, and pair creation processes 
sclf-consistently. The model spectra can fit Vela observa- 
tion well at 100 McV - 10 GeV. In their model, it was 
analytically demonstrated that an active gap, which must 
be non-vacuum, possesses a qualitatively different prop- 
erties from the vacuum solution of traditional outer-gap 
models. It was shown that the inner boundary of the outer 
gap shifts toward the stellar surface as the created current 
increases. 

Subsequently, Takata, Shibata, and Hirotani (2004) 
and Takata et al. (2006) studied the two-dimensional 
model which takes the trans-field structure into account. 
Solving the Poisson equation for the electrostatic potential 
on the two-dimensional poloidal plane, they revealed that 
the gap inner boundary is located inside of the null sur- 
face because of the pair creation within the gap. In their 
model, it is assumed that the particle motion immediately 
saturates in the balance between electric and radiation- 
reaction forces. However, there are cases in which par- 
ticles run a good fraction of the whole gap length until 
they achieve the saturation. On these grounds, we solve 
the unsaturated particle momenta and pitch angle evo- 
lution with synchrotron-curvature radiation in the two- 
dimensional geometry and we will discuss the particle mo- 
tion and the electric field structure of the outer gap in this 
paper. 

The high energy emission from millisecond pulsars 
(MSPs) is another important issue. Because of their short 
period, small period derivatives (P < 10~ 17 s s _1 ) as well 
as their evolutional history in the binary systems, MSPs 
form a distinctive group from normal pulsars. Thus, it is 
worth constructing a self-consistent electrodynamic model 
for high-energy emission from MSPs. Since the electrody- 
namic outer-gap models (Hirotani, Shibata & Takata 2004 
; Takata et al. 2006 ; Hirotani 2006a) have applied the 
method only to young and middle-aged pulsars (e.g., Crab, 
Vela, Geminga), and used whole surface cooling neutron 
star thermal emission, wc investigate the electrodynamic 



142 



Pei-Chun Hsu, Kouichi Hirotani & Hsiang Kuang Chang: Particle Acceleration in Pulsar Magnetospheres 



model of MSPs in this paper, evaluating the soft-photon 
field from the heated polar cap. 

2. Electrodynamics of Particle Accelerator 

2.1. Poisson equation 

The Poisson equation of the non-corotational potential cj> 
is given by 

w 2 cb = -Mp-p G j) , (i) 

where the Goldreich-Julian charge density is defined as 
(Goldreich & Julian 1969) 

n-B (nxBKVxB) 

PGJ = - -T, 1 a ■ ( 2 ) 



2ttc 



2.2. particle equations of motion 



Aire 



The evolutions of momentum (p = \p\) and pitch angle 
(x) of a positron and an electron are described by 



d P sc 

-jlP = qEu cos x 

at c 



dt X 



qE\ \ sinx 
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d p cos x 

dt y/m 2 e + (p/c) 2 ' 
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(5) 



where c is the speed of light. For positrons (or electrons), 
q = +e (or q = — e) is adopted, where e is the magnitude of 
the charge on an electron. Psc/c is the radiation reaction 
force due to synchro-curvature radiation. 

2.3. Two-photon pair creation rate 

The electron-positron pairs are created via photon-photon 
pair creation. The creation rate S+(t,r,p) for positrons 
and S- (t, r,p) for electrons with momentum p at position 
r and time t is given by 



S±(t,r,p) = J dfi c J de 1 J dk r] p ±g(e^,r,k) 



(0) 



where the function g represents the 7-ray distribution 
function at energy m e c 2 e 7 , momentum k, and position r. 
The pair creation redistribution function r/ p can be defined 

as 



= (1 - Pc) I 

Jet, 



de x -^-ap(e 7 ,ex,fJt c )S(p — k/2), (7) 



where fi c is the cosine of the collision angle between the 
two photons, eth = 2/[(l — Pcjt-y] is the dimcnsionless 
threshold energy for a X-ray photon having energy m e c 2 e x 
to make a pair with the 7-ray having energy m e c 2 e 7 ; a p 
refers to the pair-creation cross-section. 



2.4- "/-ray emission, propagation, and absorption 

The probability for a charged particle (with charge e) to 
emit photons via pure-curvature process is given by (Ry- 
bicki & Lightman 1979) 



1 3 hcT 3 
m e c 2 47r Rc 

F{s) = s / K 5/3 (t)dt ; 

J X 



(8) 



(9) 



(10) 



where e c is the critical energy of curvature radiation, R c 
the curvature radius of the magnetic field line, and K 5 / 3 
the modified Bcssel function of the order 5/3. 

A 7-ray photon is emitted tangential to the local mag- 
netic field line. The photon propagates along a straight 
ray path because we assume the flat space-time geome- 
try, which is a good approximation in the outer magneto- 
sphere. 

During the propagation, a 7-ray may collide with the 
thermal X-ray and be absorbed due to the pair creation 
process as described in the foregoing section. The absorp- 
tion rate is given by 



/l />oo 
dfj, c J dFr] p (r,T,fx c ) ■ g(e 7 ,r,k) 



(11) 



2.5. boundary conditions 

We consider the situation that the lower boundary of the 
gap, 8* = #™ ax , coincides with the last open field line. 6** 
is the magnetic colatitude of the point where the magnetic 
field line intersects the stellar surface. We assume that the 
upper boundary, 8* = #" lln , coincides with a particular 
magnetic field line between the last open field line and 
the magnetic axis. The thickness of the gap is defined as 



(12) 



We solve the Poisson equation from the neutron star 
surface to the outer magnctosphcre. Thus, we define that 
the inner boundary coincides with the stellar surface. As 
for the outer boundary, we solve the Poisson equation to 
a large enough distance, s = 1.4 otlc, which is located 
outside of the light cylinder. Since the structure of the 
outer-most part of the magnetosphere is highly unknown, 
we artificially set E\\ = if the distance from the rotation 
axis, U7, becomes greater than 0.90 tdlc- Under this arti- 
ficially suppressed Eu distribution in w > 0.90 tulc, we 
solve the Boltzmann equations for outward-migrating par- 
ticles and 7-rays in < s < 1.4 otlc- For inward-migrating 
particles and 7-rays, we solve only in vo < 0.9 wlc- 

To solve the Poisson equation, an elliptic-type partial 
differential equation, we impose cj> = at the inner, lower, 
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and upper boundaries, assuming that both the lower, and 
upper boundaries are grounded to the star. At the outer 
boundary, we impose d(j>/ds = —E\\ = 0. 

3. Result 

3.1. Time- dependent Simulation 

It is, unfortunately, difficult to constrain the time- 
dependent feature of 7-ray emissions observationally, be- 
cause the EGRET experiment detected one GeV photon 
every tens of thousand pulsar rotations. Taking advantage 
of this fact, all the high-energy emission theories from ro- 
tating neutron-star magnetospheres have concentrated on 
stationary analyses. This drives us to the question whether 
the particle accelerators are, indeed, stationary or not, 
which is still unsettled. 

On these grounds, we solve the Poisson equation 
and time-dependent Boltzmann equations for electrons, 
positrons and gamma-ray photons self-consistently. By 
starting from a nearly vacuum solution, we find that the 
solution evolves to a stationary state in a wide range of 
pulsar parameters, such as period, P, period derivative P, 
and magnetic inclination angle, a\. We will demonstrate 
this new result, applying the method to millisecond pulsar 
parameters in the subsections below. 

3.2. Stability of Outer- Magneto spheric Gaps 

It is conjectured by Hirotani (2006) that an outer- 
magnctosphcric gap is electrodynamically stable. In this 
subsection, we will examine this prediction by a time- 
dependent simulation. We start with the vacuum solu- 
tion of E\i(s, 0*) with an injection of electrons at 1% of 
Goldreich- Julian rate, 0.01 £IB/2tt [cm _2 s _1 ], across the 
outer boundary. 

Such injected, inward-accelerated electrons emit copi- 
ous curvature 7-rays, which efficiently collide with the 
surface X-rays to materialize as pairs in the inner mag- 
netosphere. The created pairs discharge, partially screen- 
ing the original E\ | . The discharged positrons are acceler- 
ated outwards, emitting curvature photons, some of which 
materialize in the outer magnetosphere. The discharged 
electrons, on the other hand, emit inward-directed 7-rays, 
which efficiently cascade into higher generation pairs by 
head-on collisions with the surface X-rays. 

If the created pair density becomes too large, the over- 
screened £j 1 results in a less-efficient 7-ray emission, lead- 
ing to a decrease of pairs in a dynamical time scale, ccjlc/c. 
On the other hand, if the created pair density is too small, 
the under-screened E\\ results in an increase of pairs. On 
these grounds, we can expect that an initial trial function 
evolves to the final, stationary solution after a lapse of 
several dynamical time scales. 

In figure [TJ we present the result of a stabilized solu- 
tion: The heated polar-cap temperature settles down to 



the final value kT ~ 153 eV, after about 20 dynamical 
timescales. We adopted P = 3.3 ms, P = 6.2 x 10~ 20 s s"\ 
and <Xi = 75°, with h m = 0.9. 

We also find that this solution is marginally self- 
sustained in the sense that if one of the parameter is 
changed such as P > 3.3 ms, or P < 6.2 x 10 -20 ss -1 , or 
014 < 75°, or h m < 0.9, the created current density evolves 
to a negligibly small value compared to the Goldreich- 
Julian value; in this case, the surface temperature is sim- 
ply attributed to the injected particles across the outer 
boundary and kept below ~ 50 eV. 

On the other hand, for a faster spin (i.e., P < 3.3 ms), 
or larger magnetic moment (i.e., P > 6.2 x 10~ 20 ss -1 ), 
the final stable solution, with created current comparable 
to the Goldreich- Julian current and significantly greater 
than the injected current, has a higher polar-cap temper- 
ature, kT > 153 eV. 
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Fig. 1. The heated polar-cap temperature evolution. 



3. 3. Millisecond pulsar death line 

Let us finally enlarge the electrodynamic consideration 
into the criterion for a gap to be self-sustained. If par- 
ticles are efficiently accelerated and the resultant created 
current density, j gap , is large enough, the electrons flowing 
onto the polar-cap surface will lead to a sufficient thermal 
emission from the surface. The efficient thermal photon 
flux ensures copious pair creation in the gap, and hence 
the large enough j gap self-consistently. 

However, if the external electric field, En, decreases 
by some fictitious mechanism (e.g., by artificial decrease 
of neutron-star magnetic moment, n), j gap and hence the 
polar-cap temperature, kT, decrease. The decrease of the 
surface photon flux further decreases j gap , which partially 
recovers the decrease of E\\. However, if £j| is decreased 
(e.g., by the decrease of /1) too much, even the vacuum 
solution of E11 cannot maintain the pair creation cascade 
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in the gap; in this case, j gap as well as kT tends to van- 
ish, remaining an in-active gap with the vacuum potential 
drop. 

Viewed in this light, we can define a 'death line' of the 
particle accelerator on (P, P) plane, adopting some fixed 
value of ai and h m . If pair creation is insufficient, then 
the gap upper boundary will be located near the magnetic 
pole, because the upper boundary will be formed by the 
copiously created charges that freely moves along the field 
lines and screens P|| . Thus, to find the death line, we adopt 
a large value of h m = 0.9 in this subsection. As for a*, we 
adopt two representative values, 60° and 75°. 



to a non-thermal origin due to active particle acceleration 
taking place in the magnctosphcre of a MSP, it is natural 
to interpret that the current calculation is associated with 
the Class-II MSPs. It would be possible to argue that the 
consistency of the obtained death line with the location 
of Class II MSPs on (P, P) diagram support this inter- 
pretation of their non-thermal emission as a result of the 
self-sustained magnetosphcric particle accelerator investi- 
gated in the present paper. 
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Fig. 2. P — P diagram of pulsars and millisecond death- 
line. Pulsars with period shorter than 10 ms are denoted 
by crosses. The dashed and solid lines denote the death- 
line for inclination angle 60 and 75 degree. The triangles 
denote the Class-II millisecond pulsars while the squares 
denote Class-I. 
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To determine the death line, wc investigate the 
marginal value of P below which a self-sustained gap does 
not exist for a fixed value of P. For both cases, a>i = 60° 
and 75°, we adopt three fixed values of P, 1.7 ms, 3.3 
ms, 10 ms. For 0^=60°, we find that P > 1.5 X 10" 20 
for P = 1.7 ms, P > 2.5 x 10" 19 for P = 3.3 ms, and 
P > 1.5 X 10~ 17 for P = 10 ms. For a % = 75°, on the other 
hand, P > 6.0 x 10" 21 for P = 1.7 ms, P > 6.0 x 10~ 20 
for P = 3.3 ms, and P > 4.0 x 10~ 18 for P = 10 ms. 
These conditions are plotted on (P, P) diagram in figure 
[2] with open circles. The inferred death lines from these 
six points are also shown in figure [U the dashed line de- 
notes the death line for a; = 60° while the solid line for 
a t = 75°. 

It follows from figure [5] that the three Class II X-ray 
MSPs, denoted by triangles, appear above the two death 
lines. On the contrary, Class I X-ray MSPs (squares) and 
most of other MSPs (crosses) are located below the death 
lines. Since the characteristics of Class-II members point 



